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abstract 


The  theoretical  problem  of  a  weakly  ionized,  constant 
temperature,  three  particle  plasma  in  an  externally  gen¬ 
erated  magnetic  field  is  reformulated  by  transforming  the 
set  of  14  macroscopic  plasma  equations  (continuity  and  mo¬ 
mentum  equations  for  ions  and  electrons  plus  Maxwell  >„  equa¬ 
tions)  in  14  unknowns  (ion  and  electron  number  densities 
and  velocities  plus  the  effective  electric  and  magnetic 
lelds  into  an  equivalent  set  of  four  integral  equations 
in  four  unknowns.  In  the  course  of  this  transformation, 
it  is  shown  that  plasma  behavior  can  be  interpreted  in 
terms  of  three  ambipolar  processes:  drift,  deformation 
and  diffusion.  Plasma  diffusion  is  characterized  by  twl 
diffusion  coefficients:  the  usual  Schottky  formula  ap¬ 
plying  in  the  direction  parallel  to  the  effective  magnetic 
ield  and  a  new  expression  for  the  ambipolar  transverse 
iffusion  coefficient  applying  in  directions  perpendicular 
to  the  effective  magnetic  field.  The  new  ambipolar  coef¬ 
ficient  differs  markedly  from  the  familiar  ambipolar  coef¬ 
ficient  associated  with  the  names  of  Bickerton,  Lehnert, 
Holway,  Allis  and  Buchsbaum;  and,  in  general,  it  gives  vail 
ues  for  the  transverse  diffusion  coefficient  which  are  two 
orders  of  magnitude  larger  than  those  given  by  the  latter 
It  is  concluded  that  ambipolar  diffusion  can  produce  a  transl 

verse  diffusion  coefficient  large  enough  to  account  for 
measured  diffusion  rates. 
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I.  INTRODUCTION 


A  coherent  theoretical  description  of  a  weakly  ionized 
pJasma  in  a  magnetic  field  has  been  an  elusive  goal  of  plasma 
physicists  for  many  years.  Theoretical  studies  of  three  par¬ 
ticle  plasmas  (electrons,  ions,  and  neutrals)  go  back  at  least 
to  19391,  and  these  were  preceded  by  studies  of  two  particle 
plasmas  (electrons  and  neutrals)  which  began  in  19062>3>4. 
Since  this  early  work,  scores  of  theories  have  been  proposed 
to  account  for  the  complex  and  frequently  surprising  behavior 
of  a  plasma  when  it  interacts  with  an  externally  generated  mag¬ 
netic  field.  Yet  today,  few  of  these  can  describe  even  a  lim¬ 
ited  range  of  the  observed  phenomena,  and  none  can  give  a  uni¬ 
fied  account  of  the  fundamental  processes  at  work  in  such  a 
plasma. 

Some  of  the  reasons  for  the  present  fragmented  state  of 
the  theory  can  be  found  in  the  many  restrictive  assumptions 
and  approximations  which  have  been  imposed  on  the  problem  in 
the  name  of  mathematical  expediency.  These  restrictions  have 
led  investigators  to  groups  of  isolated  and  often  contradic¬ 
tory  theoretical  results  which  fail  tc  provide  a  coherent  pic¬ 
ture  of  the  plasma.  The  best  illustration  of  this  situation 
is  the  collection  of  expressions  for  the  transverse  (to  the 
magnetic  field)  diffusion  coefficient  Dr  of  a  three  particle 
plasma.  At  least  five  expressions  for  Dfc  have  been  obtained 
under  various  assumptions: 

Dt  =  De//(1  +  Chapman  and  Cowling1  (1) 

1  For  numbered  references,  see  Sec.  V. 
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Dt  =  oikT^elBl 


Bohm  et  al.  ‘ 


Dt  =  Da /(-1  +  ^i^e5*6)  Bickerton, 6  Lehnert,7 

Allis  and  Buchsbaum, 8 
Holway9 


Dt  =  Di/(l  +  Simon, 10  Tonks11 

Dt  =  e  ^|ui/  (X  +  U-|B*  fi)  +  ^iUe/(l  +u|B‘BjJ 


Simon3 


\  / 

where  and  ne  are  the  ion  and  electron  mobilities,  d. 

and  De  the  intrinsic  diffusion  coefficients  for  ions  and  elec¬ 
trons,  t\  and  Te  the  ion  and  electron  temperatures,  and  B 
the  magnetic  flux  density.  In  addition,  a  is  a  proportionality 
constant,  k  is  Boltzmann’s  constant,  e  is  the  magnitude  of 
the  electronic  charge,  and  Da  is  the  ambipolar  diffusion  coef¬ 
ficient  of  a  plasma  in  the  absence  of  a  magnetic  field:13 
» 


Da  '  (^iDe  +  ^eDi)  /(“i  +  O 


In  this  collection,  Eq.  (1)  was  obtained  by  direct  solution  of 
Boltzmann’s  equations  for  the  ions  and  electrons.  One  might 
therefore  expect  it  to  be  more  accurate  than  the  others  which 
were  all  obtained  by  less  sophisticated  means.  in  fact,  however, 
(1)  differs  only  slightly  from  the  corresponding  expression  for 
3  tW°  Particle  Plasm.-  which  Bohm  et  al.1*  showed  was  much  too 
small  to  account  for  their  laboratory  observations.  Equation  (2) 
XB  3  semi~emPirical  relation  based  on  an  unpublished  "turbulent 
drain"  theory  of  electron  diffusion  which  Bohm  et  al.  used  to 
account  for  their  observations.  This  expression  is  unique  in 
that  it  predicts  a  fBf  variation  of  d  for  large  magnetic 
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fields.  Unfortunately ,  subsequent  experiments7 3 1 5 ’ 1 6  have 
shown  a  (B*B)-1  =  |b|  "2  variation  such  as  predicted  by  all 

classic  ( non-turbulent )  theories,  and,  thus,  the  validity  of  (2) 
is  doubtful.  Equation  (3)  is  referred  to  as  the  ambipolar  trans¬ 
verse  diffusion  coefficient  because  it  is  based  in  part  on  the 
assumption  that  ions  and  electrons  diffuse  at  the  same  rate 
across  the  magnetic  field.  Unlike  all  the  others,  this  expres¬ 
sion  reduces  to  the  correct  limit  Dfc  -»  Da  (Eq.  (6))  as  |b|*>  0. 

As  will  be  emphasized  throughout  this  paper,  the  ambipolar  as¬ 
sumption  is  generally  valid;  however,  in  obtaining  (3),  other, 
more  restrictive,  assumptions  are  imposed  on  the  particle  veloc¬ 
ities  and  internal  fields  with  the  result  that  this  expression 
(l)  gives  values  for  which  are  much  too  small  to  ac¬ 

count  for  experimental  observations7  * 1 6  except  when  u-U  B*B  -»  0. 

By  contrast,  Eqs.  (4)  and  (5)  are  in  good  agreement  with  measured 
values  of  Dt.  However,  these  expressions  were  obtained  by  aban¬ 
doning  ambipolar  diffusion  in  favor  of  a  process  in  which  electrons 
are  constrained  to  move  only  along  magnetic  field  lines  while  ions 
diffuse  across  lines  at  the  rates  given  by  (4)  and  (5).  Charge 
neutrality  is  maintained  by  currents  at  the  boundaries  which 
must  therefore  be  conductors,  either  metal  in  the  case  of  lab¬ 
oratory  plasmas  (Eq.  (4))  or  another  plasma  such  as  the  iono¬ 
sphere  in  the  case  of  a  plasma  cloud  released  in  the  upper  at¬ 
mosphere  (Eq.  (5)).  Thus,  the  derivation  of  (4)  and  (5)  de¬ 
pends  on  a  restricted  type  of  particle  motion  and  a  particular 
set  of  plasma  boundary  conditions.  According  to  this  picture, 
diffusion  should  revert  to  the  ambipolar  type  equation  (3)  when 
the  plasma  is  contained  within  insulating  boundaries.  However, 

Eq.  (3)  invariably  gives  values  for  Dt  which  are  much  smaller 
than  measured  values  even  when  experiments  are  carried  out  with 
quartz  discharge  tubes16  which  one  would  expect  to  provide  a 
working  approximation  to  a  perfect  insulator.  Thus,  the  theory 
leaves  us  in  a  dilemma  with  respect  to  Dt. 
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The  preceding  shows  that  there  is  a  clear  need  to  reform¬ 
ulate  the  problem  of  a  weakly  ionized  plasma  in  a  magnetic  field 
on  a  more  rigorous  basis  than  has  been  the  practice  in  the  past. 
The  present  paper  is  an  attempt  to  do  this  by  avoiding  res trie-’ 
tive  assumptions  and  dealing  in  a  self-consistent  manner  with 
the  macroscopic  plasma  equations  (continuity  and  momentum  equa¬ 
tions  for  ions  and  electrons  plus  Maxwell's  equations),  we  limit 
ourselves  to  an  explicit  set  of  assumptions  which  is  sufficiently 
general  to  insure  self-consistency  and  sufficiently  specific  to 
permit  the  reduction  of  the  plasma  equations  to  a  smaller  and 

more  tractable  set  of  equations.  Our  choice  for  this  set  of  as- 
sumptions  is  the  following: 

A.  inertial  forces  are  negligible  in  the  ion  and  elec¬ 
tron  momentum  equations. 

B.  Displacement  currents  are  negligible  in  Maxwell's 
equations. 

C.  Externally  maintained  electric  fields  within  the 
plasma  are  not  present. 

D.  ion  and  electron  temperatures  are  constant. 

E.  Charged  particle— neutral  particle  collision  fre¬ 
quencies  are  constant. 

F.  The  neutral  particle  velocity  is  a  known  constant. 

G.  Ion— electron  collisions  are  negligible. 

H.  Ions  are  singly  ionized  and  of  a  single  species. 

among  the  least  restrictive  of  the  conditions  usually 
imposed  in  investigations  of  this  problem  and  should  not  require 
extensive  comment  here.  Assumption  A  limits  our  considerations 
to  time  regimes  much  greater  than  the  mean  times  between  colli¬ 
sions.  B  excludes  electromagnetic  wave  phenomenon,  c  excludes 
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arbitrary  distributions  of  charge  within  the  plasma*  D  and 
E  reduce  particle  mobilities  and  diffusion  coefficients  to 
constants.  F  means  that  effects  of  ion  and  electron  motion 
on  neutral  particle  motion  will  be  ignored.  G  and  H  are 
included  as  matters  of  convenience;  both  could  be  dropped  at 
the  expense  of  encumbering  the  mathematical  details.  More  im¬ 
portant  from  the  standpoint  of  this  paper  are  the  commonly  im¬ 
posed  restrictions  which  will  not  be  used  here.  Thus,  contrary 
to  the  nearly  universal  practice  at  present,  no  restrictions 

ke  placed  on  the  motion  of  ions  and  electrons  beyond  those 
required  by  the  plasma  equations.  Hypothetical  "congruence  re¬ 
lations  8,9  will  not  be  used.  The  treatment  will  be  fully 
three  dimensional  in  the  spatial  variables;  there  will  be  no 
assumptions  of  one  or  two  dimensional  particle  distributions.17 
Steady-state  plasmas7 ’ 8  will  not  be  assumed;  partial  derivatives 
with  respect  to  time  will  be  retained  in  the  continuity  equa¬ 
tions  and,  in  general,  all  quantities  will  be  functions  of  time 
as  well  as  the  spatial  variables.  It  will  not  be  assumed  that 
the  effective  magnetic  field  within  the  plasma  is  equal  to  the 
known  external  magnetic  field  as  is  the  current  practice.  Rather, 
the  effective  field  will  be  considered  an  unknown  variable  equal 
to  some  as  yet  undetermined  function  of  the  external  field  and 
the  induced  field  generated  by  currents  within  the  plasma.  To 
do  this  in  a  self-consistent  manner  will  require  that  the  com¬ 
plete  set  of  Maxwell’s  equations  be  used.  Truncated  versions 
of  Maxwell's  equations,  such  as  those  consisting  solely  of  Pois¬ 
son's  equation,18  will  not  be  used.  Similarly,  no  preconditions 
will  be  imposed  on  the  electric  field  other  than  that  it  origi¬ 
nate  solely  within  and  from  the  plasma.  Cylindrically  symmetric 
electric  fields9  will  not  be  assumed. 

With  the  preceding  conditions  understood,  we  have  the  fol¬ 
lowing  macroscopic  equations  for  a  three  particle  weakly  ionized 
plasma  in  a  magnetic  field; 
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Continuity  equations  for  ions  and  electrons 


dm/dt  +  v-  (niVi)  =  0 

dn  /dt  +  7*  (n  V  )  =0 
^  6  0 

Momentum  equations  for  ions  and  electrons 
ni(Vi  =  -Divni  +  uini(I  +  V.XB) 

ne(ve-vn)  =  -DeVne  -  Uene(E  +  VeXB) 


Maxwell's  equations 


V  x  E  =  -dB/dt 


V  x  B  =  UQe(niVi  -  n  V  ) 


(11) 

(12) 


where  m  and  n0  are  the  particle  densities  of  ions  and  elec¬ 
trons  respectively,  V±  and  Vq  the  particle  velocities  of  ions 
and  electrons,  E  and  B  the  electric  field  and  magnetic  flux 
density,  uq  the  permeability  of  free  space,  and  V  the  neu¬ 
tral  particle  velocity.  As  previously  defined  u±  "and  y,  are 
the  ion  and  electron  mobilities,  and  D.  and  D  1  are  the^n- 
trinsic  diffusion  coefficients.  d±  and  D0  can^e  written  in 

terms  of  y,  and  y,  . 

l  e 

Di  "  kTi“i/e  -  °e  -  kTeue/e  (13) 

and  the  latter  in  terms  of  vi,  the  ion-neutral  collision  fre¬ 
quency,  and  v0,  the  electron-neutral  collision  frequency; 


T- 54/174-4-50 


-6- 


RIVERSIDE  RESEARCH  INSTITUTE 


e/m .  V . 
1  1 


M-.  =  e/m  v 
e  e  e 


(1< 


where 


k>  e,  T., 


addition,  m. 

1 


and 
and  m 


are  as  previously  defined  and,  in 


are  the  ion  and  electron  masses 


Equations  (7)  to  (12 (^comprise  a  set  of  14  equations  in 
14  unknowns  (n^.v.,  Ve,B,B) .  In  the  following  section 
thus  set  will  be  reduced  to  an  equivalent  set  of  seven  equations 

seven  unknowns.  In  the  course  of  the  reduction,  the  follow- 
ing  points  will  be  noted: 


3. 


ni  “  ne’  that  is>  strict  charge  neutrality  is  main¬ 
tained  in  the  plasma  described  by  Eqs.  (7)  to  (12) 

A  self-consistent  treatment  cannot  assume  that  the 
effective  magnetic  field  within  the  plasma  is  constant 

All  plasma  behavior  can  be  interpreted  in  terms  of 

three  ambipolar  processes:  drift,  deformation,  and 
diffusion. 


4. 


The  general  expression  for  the  transverse  ambipolar 
diffusion  coefficient  is: 


"  ^iV(l  +  ^e2B-B)  +  ^eDi/(1  +  ^i2B-B)^/(ni  +  li 


(15 


where  B  is  not  necessarily  constant. 


is 


in  Sec.  hi  a  zeroth  order  solution  for  n  *  n  «=  n 
obtained  using  initial  estimates  for  the  effective  electric 
and  magnetic  fields,  with  the  aid  of  this  solution  and  by  the 
introduction  of  a  vector  potential  the  set  of  differential 
equations  is  then  transformed  into  a  set  of  four  integral  equa¬ 
tions  in  four  unknowns  which  can  be  solved  by  successive  ap- 
proximat ions  to  obtain  improved  estimates  of  n,  E,  and  B. 
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11  *  —  REDUCED  S£T  OF  PLASMA  FIQUATIOMS 


(16) 


EquaUons  (7)  to  (12)  win  be  reduced  to  an  equivalent 
t  of  seven  equations  in  seven  unknowns  by  eliminating  one 

We  first°delnU1ty  °nS  an<J  b°th  of  the  Tnomentum  equations. 

We  first  demonstrate  that  only  one  of  the  continuity  equations 

can  be  independent  by  showing  that  „  Thus,  subtract- 

mg  (8)  from  (7),  we  have: 

^^ni“ne^^t  +  V‘ (ni^i  ”  neVe)  =  0 
But  from  (12), 

V-(n.V.  -neVe)  =  V-(V  X  B)  =  0 
Hence,  (16)  reduces  to 
^ni  “  ne)/  *  0 

Integrating  the  latter,  we  obtain: 

P(r) 


n .  -  n 
i  e 


where  p,  the  static  space  charge  density,  is  an  arbitrarv 
func  ron  of  position  only.  No„  an  arbitrary  distrl^T^ 
c  arge  could  only  be  maintained  by  externally  applied 
electric  fields.  Such  fields  are  excluded  by  Assumption  (1 

nee,  p_  °,  and  n.=ne  follow^  we  may  then  define  the 

plasma  density  n  as  follows: 


(17) 


n  =  n.  =  n 
i  e 


and  dispense  with  one  of  the  continuity  equations. 

The  preceding  shows  that  Eqs.  (7)  to  (12)  together  with 
Assumption  (c)  describe  a  plasma  in  which  strict  charge- 


(13) 
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neutrality  is  maintained.  Thus  it  n«t  - 

r  -  aS  a  rs;:r 

charge  effects  for  which  P 

"i  -  ne  “  P(r-t)  (19 

whe-e  p  ie  not  an  arbitrary  function,  cannot  be  treated  in 
self-consistent  manner  on  the  basis  of  the  plasma  equations 

:  flit  Writte"  the"  here'  In  to  treat  space 

consistently,  it  rs  necessary  to  include  the  displace¬ 
ment  term  +e  on  the  right  side  of  Eq.  (12)  j  ls  the 

presence  of  this  term  in  (12)  that  enablel  one  L  dta  h 
auxiliary  relation  the 

ni  “  ne  =  e  ve  ^  0 

(20) 

by  retracing  the  steps  from  (16)  to  (17).  m  spite  of  this 
frequent  attempts  are  made  to  treat  plasma-generated  space 

rge  effects  by  postulating  Eq.  (20)  together  with  (12)  with¬ 
out  including  the  displacement  term  in  the  latter,  Whatever 

hlh  iirerits.°f  such  a  y  „ith :::: 

Eq  \u7ZT  “  ^  '’m0n9  Si"Ce’  ^  ^  ab°TO> 

Eq.  (12)  and  the  continuity  equations  lead  inexorably  to 

era't"d  i  ^  V'E“°  in  the  absence  of  externally-ge„- 

smelted  electric  fi pihc  »  ■* 

goals  of  th  '  consistency  is  one  of  the  main 

goals  of  the  present  investigation,  we  will  not  compromise  it 

y  in  reducing  (20)  into  the  treatment.  The  following  will 
S  ow  that,  besides  being  undesirable  from  the  standpoint  of 

mi;--9-  (20)  iS  unne«ssary  for  an  adequate  theoret¬ 
ical  description  of  plasma  drift  and  diffusion. 

Following  Johnson  and  Hulburt,1’  we  can  now  eliminate  the 
momentum  equations  by  solving  (9)  and  (10)  for  the  ion  and 
electron  velocities 
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ViMXi  +  HjXiXB  +  H  2(X. -B)B-D.  2n-n.D.2!iXB 
L  1  1  1  in  i  i  n 

-  b]/b, 

ve  -  [xe  -  HeXe  X  B  +  ne2  (xe -  B)B  -  Be  a  +  ^VnjiB 


where 


Xi  =  vn  + 

Bi  =  1  +  h/b-b  , 


*e  -  v„  -  ^eE 
Be  «  1  +  |1  !B-B 


and  substituting  these  expressions  into  Eq.  (121  and  either 

of  the  continuity  equations.  In  this  way,  Eqs.  (7)  and  (12) 
become  ' 


V  XB  *  |l 


3t  +  V-^n(xi  +  M.iXxB+p12(x..B)B^-DiVn-M..DiVnxB 

VXE*v[n(VBi  -  VBe+  (^Xi/Bi  +  M.eXe/Be)  xb 

+  (<M'i2^i/Bi  •M'e%/Be)-B)Bj-(D./B.  -De/Be)Vn  ( 

-  (^iVBi  +  ^eV^7"  *B -  !Vb.  *,WBe)  (Vn-B)?] 

The  preceding  together_Jith  ( 11 !  a^a~s~et^f  lev^eq^a'tions 

in  seven  unknowns  (n,E,B)  equivalent  to  t.,e  original  plasma 
equations. 

By  solving  (26)  for  E  in  terms  of  V  x  B  and  using  the 
resulting  expression  to  replace  E  in  (25),  the  continuity 
equation  can  be  rewritten  in  a  more  convenient  symmetric  form. 
Thus,  we  find  after  a  lengthy  calculation  (see  appendix) 
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E  =  -Vn  X  B  + 


W  +  tW  n  \  ^i+^e  /  n 


+  k^f+\)5n][fe+^)  VXB 

•  XB-(-i^  ^-)(  (7  XB)  -  B)b] 

Substituting  (27)  into  (25)  and  carrying  out  another  lengthy 
calculation,  we  obtain 

3n/at  +  v.  |^nvn -DaVn+ 1x^(7  xb)  XB/tioe(|ii +(ie)J  «  o  I 

where  Da  has  been  defined  previously  (Eq.(6)).  This  form  of 
the  continuity  equation  verifies  an  observation  of  Holway20 
that  the  assumption  of  strict  "congruence"  where21 


(27) 


(28) 


n.V.  =  n  V 
ii  e  e 


(29) 


leads  to  the  conclusion  that  a  magnetic  field  has  no  effect  on 
diffusion  of  charged  particles.  For,  if  (29)  were  valid,  then 
(12)  would  reduce  to  V  X  B  *  0  and  (28)  would  become: 

Sn/St  +  v  [nvn  -  Davn]  =  0  (30) 

Which  is  just  the  diffusion  equation  for  a  field-free  plasma 
with  a  drift  velocity  equal  to  the  neutral  particle  velocity 
Vn*  Since  (30)  is  clearly  inappropriate  for  a  plasma  in  a  mag¬ 
netic  field,  it  follows  that  (29)  cannot  be  valid.  At  the  same 
time  (28)  shows  that  one  is  never  completely  justified  in  as¬ 
suming  that  the  effective  magnetic  field  is  constant  within  the 
plasma.  For  if  B  is  constant,  then  V  x  B  =  0  and  (28)  re¬ 
duces  to  (30)  which  again  is  inappropriate. 

Equation  (28)  can  be  rewritten  in  several  interesting  ways. 
Perhaps  the  most  illuminating  of  these  is  obtained  by  noting 
from  (26)  that 
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( V  X b)  XB  =  ^oe  jnCxB  -  T)VnxB  -6^(B.Vn)B  -  (b-B)Vi^J 

(31) 

where 

C  =  %  +  fE  +(p/n  +  PE)  X  B 

(32) 

i  -  Di/Bi  -  VBe 

(33) 

8  -  ^iVBi  +  ^eDe/Be 

(34) 

0  =  1/B.  -  1/B 
i  'e 

(35) 

0  =  M-. /B.  +  11  /B 
ii  e  e 

(36) 

P  ■*  ^.2/B.  -  M-  2 /b 
ii  e  e 

(37) 

and  using  this  expression  to  replace  (V  X  B)  x  B  in  (28). 
This  gives 


3n/dt  +  v-Vn+(vv)n-v[DtVn+en.ne(B-Vn)B/(n.  +  ue)l  -  0 

where 

v  -  vn  +  M-M^C  x  B/(M-i  +  M-0 )  +  V  X  (  t)B  ) 

Dt  =  (M-iDe/Be  +  ^eDi/Bi)/(^i  +  ^e) 

Decomposing  the  gradient  operator  into  components  perpendic 
ular  and  parallel  to  the  magnetic  field, 

V  «  vt  +  b(b-v)/b*b 

we  can  rewrite  the  bracketed  term  in  (38)  as  follows: 

[  ]  "  DtVtn  +  DaB(B-Vn)/B.B 
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Equation  (35)  then  becomes 


3n/dt  +  v.Vn+  (V-V)n-v.[DtVtn  +  DaB(B.Vn)/B-B]  -  0 


(43) 


,V  and  Dt  are  factions  of  the  effective  electric  end 
magnetic  fields  and  thus  subject  to  various  interpretations, 

the  °  the  COntinuit''  Ration  strongly  suggests  that 

d  ft  Tf  CaUY  threS  Pr°°esses  at  i"  the  plasma: 

drift,  deformation,  and  diffusion.  The  second  term  in  (43) 

pro  aces  an  overall  plasma  drift  velocity  of  v,  the  third 

term  accounts  for  increases  or  decreases  in  plasma  density  due 

to  contractions  (V-V  <  0)  or  dilations  (v-V  >  0)  of  the 

Plasma,  and  the  fourth  term  describes  plasma  diffusion  perpen- 

fie'lT  «  *  an<3  Parallel  (Da>  t0  thS  6ffeCtiVe  -**«ic 
Eo  S  “?  that  theSe  processes  must  be  ambipolar  because 

app  res  to  both  ions  and  electrons,  (it  does  not  mat¬ 
ter  whether  one  starts  with  the  electron  continuity  equation 
(8  or  the  ion  continuity  equation  (7)  as  was  the  choice  here- 
both  of  these  equations  lead  to  (43)  using  the  method  employed 
in  this  section.)  We  note  further  that  D  (Eq.  (40))  ap- 

proaches  the  correct  limit-  n  r 

correct  limit  Da  as  B  goes  to  zero  and  that 

tor  nonzero  B 
D.  s*  D./b. 

1  1  (44) 

provided  the  usual  ordering  ue  »  m-j.  holds.  This  approxima¬ 
tion  13  ldentl«l  the  expression  obtaiued  by  Simon  (Eq.  (14)) 
but  it  is  not  restricted  by  the  latter's  assumptions,  it  shows  ’ 
that  special  restrictions  need  not  be  imposed  on  the  ion  and 
e  ectron  motions  and  that  ambipolar  diffusion  can  indeed  pro- 
uce  a  transverse  diffusion  coefficient  large  enough  to  account 
or  measured  diffusion  rates.  Finally,  it  should  be  emphasised 
that  no  restrictions  have  been  placed  on  8  in  obtaining  (40) 
an  (44).  Thus,  these  equations  can  be  applied  on  a  point-by- 
point  basis  in  the  general  case  when  8  is  not  constant. 


T- 54/174-4-50 


-13- 


RIVERSIDE  RESEARCH  INSTITUTE 


111  •  TRANSFORMATION  OF  THE  PT.ASMA  EQtlATTAMg 
INTO  A  SET  OF  INTEGRAL  EQUATIONS 


Equations  (39),  (40),  and  (43)  provide  a  basic  descrip. 
tion  of  the  plasma  in  terms  of  ambipolar  drift,  deformation 
and  diffusion;  and,  as  such,  these  equations  contain  the  most 
important  results  of  the  present  paper,  obtaining  an  actual 
solution  to  (43)  for  the  plasma-density  distribution  in  con¬ 
junction  with  Maxwell's  equations  ((11)  and  (26))  remains  a 
formidable  problem  which  will  not  be  attempted  here,  m  this 
section  we  will  limit  ourselves  to  describing  one  approach  to 
the  problem  which  is  currently  under  investigation. 

This  approach  is  based  on  transforming  Eqs.  (11),  (26), 
and  (43)  into  an  equivalent  set  of  integral  equations  which, 
in  principle  at  least,  is  amenable  to  solution  by  successive 
approximations.  To  do  this,  we  note  first  that  if  the  electric 
and  magnetic  fields  were  constant  say  E  and  B  ,  then 
Eq.  (43)  would  reduce  to:  °  ° 

3no/3t  +  •  Vno  -  DtoVto*„o  -  DaBo  •  V(B0  .  Vno)/r  .  Bq  .  0  (45) 

Where  VQ  and  _Dto  are  given  by  ( 39)  and  (40)  with  E  and 
B  repiacing  E  and  B.  Equations  of  this  form  can  be  solved 
exactly.  The  solution  to  the  initial  value  problem  for  N 

particles  distributed  in  the  form  of  a  spherical  Gaussian  with 
radius  Rq  at  t  =  o  is: 


n  (r,t)  = 


7T3//?LT2 


exp  - 


O  o 


(46) 
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where  we  have  decomposed  the  position  vector  r  and  drift 

velocity  Vq  into  components  perpendicular  and  parallel  to 

B  : 
o 

7  =  7to  +  *o<V7)/V*o 


V 


=  Vto  +  Bn(Bo‘Vo)/B.-B, 


and  we  have  defined  T  and  L  as  follows 


T2  «  4d.  t  +  R  2 
to  o 

L2  =  4D  t  +  R  2 
a  o 

Equation  (46)  describes  a  cylindrically  symmetric  particle 
distribution  moving  with  a  constant  drift  velocity  V  and 
simultaneously  diffusing  perpendicular  and  parallel  to  B 
with  characteristic  diffusion  coefficients  D  and  D  ° 
respectively.  The  contours  of  constant  density°are  ellip¬ 
soids  of  revolution  with  major  axes  parallel  to  B  .  Now, 
in  general,  E  and  B  are  not  constant,  and  thus°we  cannot 
expect  (46)  alone  to  give  an  adequate  description  of  the  man¬ 
ner  in  which  an  initially  spherical  plasma  distribution  evolves 
in  a  magnetic  field.  Indeed,  observations  of  barium  plasma 
clouds  in  the  upper  atmosphere22  show  that  the  ultimate  dispo¬ 
sition  of  an  initially  spherical  plasma  distribution  is  much 
more  complicated  than  (46),  involving  the  development  of  mag¬ 
netic  field-aligned  striations  among  other  things.  Neverthe¬ 
less,  (46)  can  provide  an  acceptable  description  of  the  early 


(47) 

(48) 

(49) 

(50) 


time  behavior  of  such  a  plasma,  if 


Bq  and  Eq  are  properly 


chosen.  For  example,  if  one  sets  Bq  equal  to  the  ambient 
geomagnetic  field  and  tabes  Eq  to  be  the  first  term  in  (27) 


E  =  -  VX  B 
o  no 


(51) 
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then,  substituting  (51)  into  (39),  we  find  that  the  drift  ve- 
locity  reduces  to: 

V  B  v 

°  r  (52) 

and,  in  this  case,  (46)  does  provide  a  reasonable  approximation 
to  what  usually  happens  during  the  first  two  or  three  mintues 
following  a  barium  plasma  release  at  heights  bstween  100  and 
300  km.  That  is,  the  plasma  cloud  initially  tends  to  drift 
With  the  velocity  of  the  ambient  neutral  particles  while  dif¬ 
fusing  parallel  and  perpendicular  to  the  geomagnetic  field  at 
rates  controlled  by  Da  and  Dto>  respectively,  the  latter 
being  determined  from  (40)  with  B  set  equal  to  the  geomag¬ 
netic  field  at  the  height  of  release.  We  may  interpret  (51) 
and  (52)  in  reverse  order  by  saying  that  initially  the  neutral 
particles  in  colliding  with  the  charged  particles  tend  to  drag 
the  plasma  at  the  neutral  particle  velocity  v  thereby  in¬ 
ducing  an  electric  field  given  by  (51).  Eventually,  other 
fields  tending  to  oppose  the  plasma  motion  will  be  induced  in 
the  cloud  resulting  in  a  smaller  drift  velocity  and  a  much  more 
complicated  plasma  distribution,  in  spite  of  this,  the  success 
o(6  in  describing  the  early  time  behavior  of  barium  plasma 
clouds  suggests  that  this  expression  is  a  valid  teroth  order 
solution  which  could  be  used  as  the  first  step  in  a  successive 
approximation  solution  to  Eqs.  (11),  (26)  and  (43).  To  facil¬ 
itate  such  a  solution,  we  can  add  v -Vn-D  V  *n-D  B  .V(B  -Vn)/ 
B0.B0  to  both  sides  of  Eg.  (43)  and  rewrite  («)  as  ?oUows? 


3n/dt  +  VQ.Vn  -  DtoVto*n  -  Vo'’<  V^/V^o  m  q[„] 


where  the  operator  Q(-l  is  defined  by 
*  V*.VM  +  (v*v# )  [•  ]  +  V*F[-] 


(53) 


(54) 
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V  =  V  -  V 
o 


(55) 


F[0-»tvt[.]-»tovto[.]  +  Da(ili2Li> 

\  B  *  B  B o  •  B  J 

ZZ  T  T  V  Green'S  the°rem'!3  (53>  »“  »•  transformed 
into  an  integral  equation  for  n: 

t  00 

"(t.t)  -  no(r,t)  +  i  Jdt>  J{Jd7'3Q[n(7-,t')]G(?-7.,t-t') 

O  —oo 

where  „0(F,t)  is  given  by  ( 46) ,  G(r, t)  is  defined  by: 


(57) 


G(r,t) 


■ - ail! -  exp  f.  (rt°-Vtn! 

,(4Dtot)(4Dat),/»  [  4D.  t 


7r3/2(4Dtot)(4Dat 

!V(?-lt))‘| 

(V5o)4V  J 


(58) 


and  in  the  latter,  H(t)  is  the  step  function.  The  form  of 

iith  "n  (Fat  7  TTV  SOlUti°n  by  aUCCeaaiva  approximations 
0(  .  )  as  the  first  approximation  to  n.  in  such  a 

process,  it  is  necessary  to  calculate  the  effective  electric 
t'on^Tboth  fleldS  a‘_eaCh  iteration,  since  Q[.]  is  a  func- 
weu.  ecu  l  easi1^  ^°ne  when  Max- 

the  v  !  nS  are  W-itten  in  integral  form‘  introducing 

the  vector  potential  A  in  the  usual  way 


B  «=  V  X  A  +  B 


(59) 


and  employing  the  coulomb  gauge  V-A  .  0  (since  „.-„)  w„ 
find  that  Maxwell's  equations  reduce  to:  1  * 
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E  =  -  dA/dt  +  E 

O 

and 


( 60 ) 


V2A  -  -  M-QJ[n,A] 

where  M-QJ[n,A]  is  given  by  the  right  side  of  (26)  with  (59) 
and  (60)  replacing  B  and  E.  Equation  (61)  may  be  consid¬ 
ered  a  vector  Poisson  equation  with  a  formal  solution24 


A(r,t) 


J[n(r',t),A(r ',t)  j 

|r -r  '  | 


(62) 


Since  j  depends  on  A,  as  well  as  n,  Eq.  (62)  is  actually 
a  vector  integral  equation  (equivalent  to  three  scalar  integral 
equations)  which  together  with  (57)  comprises  a  coupled  set  of 
nonlinear  integral  equations  in  the  unknowns  n  and  A.  This 
set  is  exactly  equivalent  to  the  original  set  of  macroscopic 
plasma  equations,  and  its  solutions  can  be  used  with  Eqs.  (18), 
(21),  (22^,^(59)  and  (60)  to  recover  complete  solutions 

(ni,ne,E,B,Vi» ve)  to  Eqs*  (7)  to  (12)»  We  have  thus  succeeded 
in  replacing  our  original  set  of  14  nonlinear  differential 
equations  in  14  unknowns  with  an  equivalent  set  of  four  non¬ 
linear  integral  equations  in  four  unknowns. 

The  form  of  Eqs.  (57)  and  (62)  suggests  a  solution  by 
successive  approximations  with  n  =  nQ  and  A  =  0  as  initial 
estimates.  The  success  of  this  approach  will  depend  on  the 
proper  choice  of  Bq  and  Eq  among  other  things,  in  the  case 
of  barium  plasmas  in  the  upper  atmosphere,  preliminary  calcu- 
us*ng  the  aml>ient  geomagnetic  field  for  B  and 
~Vn * Bo  for  Eo  are  enc°uraging. 
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IV.  DISCUSSION 


In  the  preceding  sections  we  have  reformulated  the  theo¬ 
retical  problem  of  a  three  particle  plasma  in  a  magnetic  field 
by  reducing  the  set  of  macroscopic  plasma  equations  to  a  smaller 
equivalent  set  of  equations,  one  of  these,  the  continuity 
equation  (43),  consists  of  terms  describing  three  ambipolar 
processes:  drift,  deformation,  and  diffusion.  From  this  equa¬ 
tion,  two  characteristic  diffusion  coefficients  emerge:  the 
Schottky  formula  (6)  applying  for  diffusion  parallel  to  the 
magnetic  field  and  Eq.  (40)  applying  for  diffusion  perpendic¬ 
ular  to  the  magnetic  field.  Equation  (40)  differs  markedly 
from  the  currently  accepted  expression  for  the  transverse  ambi- 
po.ar  diffusion  coefficient  (3),  and,  in  general, (40)  will  give 
a  much  larger  value  for  D^ .  To  show  this,  we  form  the  ratio 
R  of  (40)  to  (3)  as  follows: 


~  (1  +  ^i^eB-B)/2(l  +  M-,2B.B) 


In  most  ecses  of  interest,  \l  al  b-B  »  1  and  [i.2B.B  »  1- 
hence,  1  ' 

R  a:  |jL  /2  M- •  . 
e  i 

Or,  using  (14), 

R  ae  m .  v  .  /2m  v 
i  i  e  e* 

For  an  02+  plasma,’  mi/me  =  6  x  10“  and  v . /v  =0.5  10'2 
Hence,  R=150.  That  is,  for  a  typical  plasma1,  the  transverse 
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diffusion  coefficient  computed  from  (40)  will  be  larger  by  two 
orders  of  magnitude  than  that  computed  from  (3).  since  the 
larger  diffusion  coefficient  is  in  substantial  agreement  with 
measured  values  of  Dt  and  since  this  coefficient  is  inde¬ 
pendent  of  any  boundary  conditions,  it  follows  that  Eg.  (40)  is 
a  means  of  escaping  the  dilemma  described  in  the  introduction. 

Equation  (40)  characterizes  transverse  diffusion  in  a 
plasma  wherever  charge  neutrality  is  maintained  as  prescribed 
by  Eqs.  (7)  to  (12).  Indeed,  one  may  interpret  (40)  as  charac¬ 
terizing  the  diffusive  properties  of  a  magnetized  plasma  con- 
sisting  of  charged  particles  which  are  constrained  to  move  in 
such  a  way  that  the  average  net  charge  per  unit  volume  is  al¬ 
ways  zero,  one  would  not  necessarily  expect  this  expression 
to  apply  in  regions  where  significant  deviations  from  charge 
neutrality  occur,  such  as  in  the  positive  ion  sheath  which 
forms  when  a  conductor  is  immersed  in  the  plasma.  As  was 
pointed  out  in  sec.  II,  this  more  general  case  can  be  treated 
in  a  self-consistent  manner  only  by  including  the  displacement 
term  in  Maxwell's  equations.  Nevertheless,  in  most  cases  of 
interest,  charge-neutrality  is  maintained  throughout  the  greater 
part  of  the  volume  occupied  by  the  plasma,  and,  in  these  cases, 

(  )  will  provide  a  means  of  computing  D^. 
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VI.  APPENDIX 

DERIVATION  OF  E QUA T TDM  (21' 


lows : 


With  the  aid  of  (23),  Eq.  (26)  can  be 


rewritten  as  fol- 


VXB  -  M.oe[n(^Vn  +  ^E+(^Vn  +  PE)  XB+((PVn+CE).B)B 
-  T]Vn  -  evn  X  B  -  a(Vn-B)Bl 


where : 

0  =  1/B.  -  1/b 
i  e 


*  “  V'Bi  +  \\ 
e  =  ^Vb.  -  neVBe 

C  -  +  M.e3/Be 


T1  =  d./b.  -  d  /b 

ii  e  e 


e  =  M.iD./Bi  +  M-eDe/Be 

°  =  W8!  -  W*. 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 


(A6) 


(A7) 


Rearranging  (Al),  we  have: 


0 E  +  PE  X  b  +  £ (e*b)r  -  Ti 


(A9) 
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where: 


3-  1*1 


+  n  a  +  e  aua  + 


-  «Vn  -  nn  X  8  -  P(Vn.B)B 


Hence, 


I'B  =  P.  -  pE  X  B  -  C(E.B)B 

,  Forming  the  dot  product  of  (All)  with  b,  we  obtain, 
^E-B  *=  fi.B  -  C(E.B)(B.B) 

Hence,  solving  for  e\ 3" 

E-B  *=  (fi-B)/(^+  £b*B) 

Similarly,  taking  the  cross  product  of  (All)  with  B, 

^XB.flxB-p(i.B)xB 

or 

K'EXB-SxB  -  g(B(E.B)  -  E(B.B)  ) 

Substituting  (A13)  into  (A15),  we  obtain: 


E  X  B  m  I  {2  X  B 


■•(■tSh)- *'•'>)]/- 


And,  substituting  (A13)  and  (A16)  into  (All),  „e  obtair 


t'E  «  ft  -  p  I  Q  x  B  -  P 


/jlisla. .  E(- -  \1  1IILSlI 
Vk'+cb.b  )\ 


(A10) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 


T- 54/174-4- 50 


-24- 


RIVERSIDE  RESEARCH  INSTITUTE 


Collecting  terms 

_-_BjBe»VXB  B.BeP(VxB)xB  (  V  X  B)  .B)B 

(•1i  +  %)2Vn  "  K  +  M-elV" 

+  B.B.(^e-eB)vnxB 

- - +  — - _ _ _ _ _ 


+  — ; - r, — + 

(^i  +  M-e)  2n  (ix.  +u)»n 


.  'W^tB-B))^  B ^(4 

\  U'i+>xe)J  (ixi  +  ne)  /  „ 


,0(0  +  3B  •  B )  Vn 


BiBe(0‘  -  30)Vn  XB  m-.m- 

-  — -  _ _ _  .  1  i 


+  >xe) 


(M-i+IXe)2  +(,.+,e) 


[l  ,  |JL  / 

-^y(0  +  @(B.B))(Vn.B)B  (A21) 


NOW 

BiBe(^TH-P0B,B) 


+0 


(|Xi  +  |Xe) 


PL.D. 

l  l 

LBi2 


M-  D  p.  D. 

e  e  +  e  i  _  __ 

V  BeB.  B-B 


itif e 

eBi 


+  7f(Bi-1)-7T(Be-1)  -^i(Be-l)+^ 

Bi  B3  BiBe  BiBe 


(Bi  -  1) 


“M-.D  +11.0.  -|i  D 
-  e  1 _ l  e  j.  i  e  e 


And 

B1Be(^e  -Pn)  B 

( M-i  +  M-e )  J  (^i  +  ^e) 


(M-i+  ^e>a 


D.  -D 
i  e 

^  +  H 

i  e 


(A22) 


h?e  rUe^iDi+^eDe  +  ^i2De  +  ^e2Dil 
+  "e>2  L'  V.  J 


¥e  +  ^ 

Hi+He 


(A23) 
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Also 

Bj Be(^  -P0) 
(  +  M-e )  2 


And 

BjB^(0+  0b.B) 
(^i+^e)2 

Furthermore 

BiBe(»a  +  P0) 

K+^J2 


^i^eCn  +  aCB.B)) 


^i+^e 


eDe 


Ox,  +  ^)>  [Vl*Di  *  •V.’D.  +  li.V,  -  \ 

’  +  ^e^^  +  aB-B)] 

^i/e  f  D.  D  /m-.2D  ll  2n  \  n 

“i+U  1  ‘  De  -  B~  +  r  '  (t^  -  -Vs)  b.b] 


D. 


Ui+ligL  i 

LI  - 

e 

B.  + 
1 

Be 

X  X 

B.  '  +  ' 

1 

.  BiBe» 

/_L 

1 

^e2\  -  _\ 

(Hi  +  ^l2 

(Bi 

"  Be 

ll7  ' 

v BB) 

_  BiV 

/-i. 

1 

.  Bi‘1 

Be-n 

vBi 

"  Be 

+  B. 

1 

B  /  “ 

e  ' 

(A24) 


(A25) 


-fA  _  .  /V  V 

A  +  V’  VB,*  Bfi2  B,Bfi 

^4-^  +  ^l  +  AL\ 

Bi  Be  BiBe  BiBe/ 

-i^g_/4i2  +  ^2  +  2^n  \ 

[ M^i  +  M-  )  2  \  B.B  / 

i  e 


(A26) 
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With  the  aid  of  (A22),  (A23),  (A24),  (A25)  and  (A26),  (A21) 
becomes: 


which  is  Eq.  (27). 
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